Boson-like quantum dynamics of association in ultracold Fermi gases 
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We study the collective association dynamics of a cold Fermi gas of 2N atoms in M atomic 
modes into a single molecular bosonic mode. The many-body fermionic problem for 2 M amplitudes 
is effectively reduced to a dynamical system of min{iV, M} + 1 amplitudes, making the solution 
no more complex than the solution of a two-mode Bose-Einstein condensate and allowing realistic 
calculations with up to 10 4 particles. The many-body dynamics is shown to be formally similar to 
the dynamics of the bosonic system under the mapping of boson particles to fermion holes, producing 
collective enhancement effects due to many-particle constructive interference. Dissociation rates are 
shown to enhance as the number of particles whereas association rates are enhanced as the number 
of holes, leading to boson-like collective behavior. 
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The role of Bose stimulation in the dynamics of coupled 
atomic and molecular Bose-Einstein condensates (BECs) 
has been extensively studied since the first theoretical 
proposals for the conversion of atomic condensates into 
molecular BECs 0, H H H @. It has been shown 
that the entire condensate will undergo large amplitude 
coherent oscillations between atoms and molecules. The 
Bose-enhanced oscillation frequency is predicted in mean- 
field theory to scale as y/~Ng where g is the single-particle 
atom-molecule coupling frequency and TV is the total 
number of condensate particles. Quantum- field effects 
slightly modify the effective frequency to (V~N / In N)g 
and introduce collapse and revival of the coherent oscil- 
lations due to interparticle entanglement 0, 0] • Yet de- 
spite these modifications the collective behaviour remains 
significantly different than the single pair dynamics, in- 
dicating the dramatic effect of quantum statistics on the 
many-body dynamics. 

Experimentally, the effort towards a molecular BEC 
and the realization of the theoretical predictions, was 
made using the techniques of Raman photoassociation 
and Feshbach resonance magnetoassociation pj.l9l.licj. 
Unexpectedly, the first molecular condensates were pro- 
duced using the Feshbach method from nearly degenerate 
gases of fermion atoms 0, 0, ^| which are more stable 
against vibrational quenching than their bose-atom coun- 
terparts. Molecular BECs made of boson constituent 
atoms came a close second 0, . 

The possibility of coupling a nearly degenerate atomic 
Fermi gas with a molecular BEC raises interesting ques- 
tions about the nature of the ensuing collective dynamics. 
At first site it appears that the dynamics of associating 
fermions should be drastically different than the associa- 
tion dynamics of bosons, since the former are subject to 
Pauli blocking as compared to the Bose stimulation af- 
fecting the time evolution of bosons. However, it has been 
previously noted that collective effects in fermion sys- 
tems can mimik bosonic stimulation in four-wave mixing 
\Hi Ell , where the boson- like behaviour is attained from 
various pathways adding up constructively [l6j . More- 
over, few-particle numerical results for the association 
of a Fcrmi-Bosc mixture of atoms Il8|, indicate some 



similarity with the purely bosonic case, in that collec- 
tive Rabi-like oscillations and Rapid adiabatic passage 
are observed in both cases. Better understanding of this 
boson-like behaviour, allowing for the simplification of 
calculations and illuminating its origin is highly desir- 
able. 

In this work we study the association dynamics of an 
atomic fermion gas into a degenerate molecular bose gas. 
We show that the many-fermion problem which is appar- 
ently much more obstrusc than the corresponding two- 
mode bosonic problem (or than the Fcrmi-Bosc mixture 
case), can be effectively reduced to an (N + 1) x (N + 1) 
system of dynamical equations, where N is the total num- 
ber of pairs (i.e. the sum of the molecule number and half 
the atom number) . The complexity of the system is thus 
identical to that of the dynamical equations obtained for 
two-mode bose association in its number-state represen- 
tation 0, Q , allowing for numerical calculations with up 
to N = I0 4 particles with current computation power, 
as opposed to characteristically N = 10 particles in Ref. 







Furthermore, we show formally that in the limit of 
large N and M where the number of fermion modes M 
is equal to the number of particles N, the structure of the 
fermionic system is a 'mirror image' of the bosonic sys- 
tem, producing precisely the same dynamics for fermion 
dissociation as was attained for boson association. Sim- 
ilarly the dynamics of association into fermions repro- 
duces the known results for boson dissociation, including 
its modulational instability It is evident from our 
model that the enhancement effects come from adding 
up the various pathways connecting states with N — n 
molecules and n dissociated pairs with states containing 
one more (or one less) molecule and one less (or one more) 
dissociated pair. Thus, the origin of the nonlinear collec- 
tive behaviour is found to be identical to the four-wave 
mixing case |lfjj . 

We consider the association of a Fermi gas of atoms 
into Bose molecules. The interaction representation 
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FIG. 1: Collective dynamics of association of an atomic 
Fermion gas into a molecular BEC (a) and dissociation of a 
molecular BEC into a degenerate Fermi gas of atoms (b). Par- 
ticle pair numbers are fixed to N = 5 (dashed line), N — 50 
(dotted line), and TV = 500 (dash-dotted line) pairs. The 
number of fermion modes is set equal to the number of parti- 
cles. 
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where e k = h 2 k 2 /2m is the kinetic enery for an atom 
with mass ra, £ k is the molecular energy containing both 
kinetic and binding conributions, c k (T are the annihila- 
tion operators for the atoms, obeying the usual Fermi 
anticommutation relation, 6 k are the molecular annihila- 
tion operators, commuting with c k .o- and obeying bosonic 
commutation relations, and g is the strength of the atom- 
molecule coupling. Since the pairing instability occurs 
predominantly for colliding atom pairs whose center of 
mass momentum is zero as confirmed by the for- 
mation of a molecular BEC in the experiments of Refs. 
[ill ll^. Jli l j w e use the single molecular mode approxi- 
mation [2fJ,|21|, assuming that only singlet pairs of oppo- 
site momenta are associated, resulting in molecules in the 
lowest energy molecular mode. Moreover, since the Fermi 
energy is characteristically very small compared with the 
atom-molecule coupling frequency, we can safely ignore 
the variation in e k . This approximation is tantamount 
to assuming that the motional timescale of the atoms is 
much longer than the timescale of atom-molecule conver- 
sion [lsj and is verified by the nearly complete conversion 
of Fermi atoms to Bose molecules in Feshbach sweep ex- 
periments 0, 0, E| (we note paranthetically that the 
variation in e k can be easily incorporated into our model 



by proper renormalization of the amplitudes), 
suiting rotating wave Hamiltonian is 
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where A = 2ek F — £q- In what follows we shall assume 
that A = 0, i.e. we consider the case of resonant coupling 
between atoms and molecules, retaining only the inter- 
action Hamiltonian. Since the Hamiltonian J5J) conserves 
the number of atoms, i.e. 2n + k n CTjk = 2N where n 
is the numeber of molecules, nk,<r is the number of atoms 
in state k, a and N is the (conserved) total number of 
pairs, we can fix N and expand the quantum state of the 
system in pair number states: 
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where k denotes some labeling of the M available momen- 
tum states, rife = 0, 1 are fermionic occupation numbers 
for the free atom states (where = 1 denotes a singlet 
pair of particles occupying the momentum states k, — k) 
and the summation over {n^} for any given n is over 
all nonordered combinations of n atom pairs among M 
states, so that X)&=i n k — n - The effective phase-space 
for the problem can thus be classified into min{iV, M} + 1 
manifolds characterized by the number of dissociated 
pairs n = 0,1,... , min{iV, M}. For each n, the size of 
the corresponding manifold is the number of orderings 
of n pairs in M states given by the binomial coefficient 
(n 1 ). It is immediately seen that the interaction Hamil- 
tonian couples each state in the manifold with n disso- 
ciated pairs to M — n states having n + 1 dissociated 
pairs and to n states with n — 1 free atom pairs. Substi- 
tuting the expansion (J3J) into the Schrodingcr equation 
ihdt\ip) = H\ip), with the Hamiltonian JTJ, we obtain: 

iC'jv-n,m,...,n M = g (\/N - (n - 1) (4) 

M 

E ^ x < n k^N-(n-l),m,...,n k -l,...,nM 



+VN - n 

m \ 

x E ^Q,n fc CW-(n+l),m,...,nfc+l,...,nM J 
k=l / 

where n = 0, . . . , min{.ZV, M }. These equations of mo- 
tion can be put into the matrix form: 

iC n = g (y/N - n + lflr'C^, + x/F^ J D»+ 1 C n+1 ) , 

(5) 

where C n (t) is a column vector of all the (^) amplitudes 
corresponding to the possible arrangements of the n atom 
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pairs among the M available fermion modes. The matri- 
ces Dj are (f 1 ) x (j) dimensional and contain only unit 
and zero elements. Their explicit form is determined by 
the ordering of elements in the vectors Cj and Cj. 

The system (JHJl contains equations for 2 M amplitudes, 
making effective dynamical calculations restricted to lim- 
ited numbers of particle (characteristically no more than 
20 particles with current computation power). How- 
ever, it is clear that the high symmetry of the system 
should produce conservation laws that would allow its 
simplification. In fact, since all the states corresponding 
to various atom ordering within each (?f) manifold of 
N — n molecules are equivalent (their coupling strength 
to states with one more or one less molecule are equal 
and each is coupled to the same number of states above 
and below it), there are — 1 constraints in each such 
manifold, leaving effectively only min{iV, M} + 1 free am- 
plitudes. Consequently, an enormous reduction in the di- 
mensionality of the system is attained by multiplication 
of Eq. JSJ) by a row vector u n> M of (^ f ) unit elements 
(thus summing over all amplitudes in the manifold of n 
dissociated pairs) and use of the identities: 



um/A" 1 = (M — n+ l)u„_i >M . 
The resulting equations of motion read 

ia n = g N — n + ly/n(M — n + l)a n - 
+VN - ny/(M -n){n + l)o<„+i) 

where 



n+1 _ 



(n + l)Un+l,M , 
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(9) 

is the sum over all (^) amplitudes with N — n molecules 
and n free atom pairs, corresponding to the various order- 
ings of the atoms, normalized by the number of permu- 
tations of n atoms in M states. The average number of 
molecules, provided that all amplitudes in the n-th man- 
ifold are equivalent (i.e. under the restriction that initial 
conditions do not break this equivalence), is calculated 
according to 

min{JV,M} 

(b$bo) = E (tf-^E \C N -n, ni ...,n M \ 2 
m=0 {n k } 

min{N,M} 

= £ (^-«)l««| 2 (10) 
ra=0 

The size of the system (jHJ) is thus reduced to 
min{iV, M} + 1, exactly as would be obtained for the 
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FIG. 2: Collective dynamics of association (a) and dissocia- 
tion (b) in a two-mode atom- molecule BEC. Particle numbers 
are fixed to N = 5 (dashed line), N = 50 (dotted line), and 
N = 500 (dash-dotted line) pairs. The number of fermion 
modes is set equal to the number of particles. 



equivalent two mode bosonic system, allowing numerical 
results with up to 10 4 particles. The physical interpreta- 
tion of Eqs. (JHJ is quite striking. It shows that in addition 
to the usual bosonic factors of \/N — n + 1 for associa- 
tion and y/N — n for dissociation, originating from the 
molecular bosonic field operator bo, the association rate 
depends on the number of pairs in the initial state (a 
trivial factor originating from having n pairs which can 
associate) but is also enhanced by the number of holes in 
the target state. Similarly, the rate of dissociation (i.e. 
the annihilation of a hole) scales with the number of holes 
in the initial state but also with the number of particles 
(or more precisely, particle-pairs) in the target state. 

The resulting dynamics should thus be surprisingly 
similar to the dynamics of an atom-molecule condensate. 
It is illuminating to compare Eq. |jSJ to the (N + 1) x 
(N + 1) number-state representation of the Schrodinger 
equation for a two- mode atom- molecule BEC 0, 



iPn = g (VN -n+ ly/2n(2n - l)/3„_i 

+VA^V(2" + 2)(2n + l)/3„+i) , (11) 



. . . , TV and the average molecule number is 
Q (N — n)\/3 n \ 2 . Examination of Eq. lfTT|) 



where n = 0, 
m = En= 

shows that both the association rate and the dissocia 
tion rate are enhanced as the number of atoms in the 
target state. Moreover, it is evident that up to factors 
of two, the main difference between the fermion equation 
|jHJ and the boson equation (|11J) is in the enhancement 
factors where particle number terms n in the boson case 
are replaced by hole number terms M — n for fermions. 
It is therefore clearly evident that when N = M (i.e. the 
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number of pair states is equal to the number of pairs) the 
dynamics of fcrmion association (dissociation) should be 
qualitatively very similar to the dynamics of boson dis- 
sociation (association), up to a factor of two in the per- 
tinent rates. 

This observation is confirmed by numerical integration 
of Eq.® and Ea. l|ll|l for fcrmions and bosons respec- 
tively. In Fig. 1 and Fig. 2 we plot the expectation value 
of the pair-number difference w = 2(6q6o) /N—l as a func- 
tion of the rescaled time r = \ // 2Ngt. Initial conditions 
had all free fermion atoms in (a) or all molecules in (b). 
The fermion results in Fig. 1 mirror the boson results 
of Fig. 2 (identical to Refs. 0,0) in that the dynam- 
ics of fermion association closely resembles the dynamics 
of boson dissociation and vice versa. The similarities 
include the y/N /InN scaling of the atom-molecule oscil- 
lation frequency [{| and the dynamical instability of the 
fermion atoms to molecule formation, noted previously 
in mean- field studies ji^, ■ This instability turns out 
to be the mirror image of the familiar modulational in- 
stability of the molecular mode in the two-mode BEC 
dynamics [f| . As expected, a factor of two exists in the 
characteristic frequencies. 

The mirror symmetry between fermion and boson as- 
sociation dynamics in the limit of large N = M is in fact 
exact. To show this we note that for large N the boson 
equations of motion for the amplitudes Pn-th may 
be written as 

N - n = 2g(Vn+l(iV-n)/3jv-„-i 

+V^{N - n + l)/3iv-„+i) , (12) 

where we have neglected some terms of order 1 <C 2N. 
Mapping /3j\r— n ~ * Pn we obtain that 

i$ n = 2g((N-n + l)VEp n - 1 

+(N - n)V^TTp n+ i) , (13) 

which is identical up to a factor of two, to the fermion 
equation (jHJ when N = M. This factor of two emanates 
from having a single atomic mode with 2n particles in 
the bosonic problem, making it the analog of a degenerate 
optical parametric amplifier (OPA), as opposed to having 
n particles each of spin f and spin J, atoms in the fermion 
problem. The correspondence is thus even better with 
the nondegenerate OPA-like association of a two species 
BEC 

H = g (aia 2 tf + ?>4 a i) > ( i4 ) 

where a\,a 2 and b are all boson annihilation opera- 
tors. The number-state representation of the Schrodingcr 
equation with the Hamiltonian (|14|l maps exactly to Eq. 
(|13|l under (3n-u — * ftm without the multiplying factor 
of two and the restriction of large N. The resulting dy- 
namics of this boson problem at any N is thus just the 
mirror image of Fig. 1. 



The dynamics of the fcrmion amplitude corresponding 
to n atom pairs (and M — n hole pairs) thus maps pre- 
cisely for N = M, to the exact quantum dynamics of 
the two-species boson amplitude corresponding to N — n 
atom pairs. At the limit of large N, the two-species 
boson dynamics reproduces the single atomic mode bo- 
son dynamics with a factor of two in the characteris- 
tic timescales. One can thus solve the fully quantum 
fermion problem by simply solving the nondegenerate 
(two atomic species) two-mode BEC problem, and map- 
ping boson amplitudes with n particles to fcrmion ampli- 
tudes with n holes, or at the limit of large N, solve the de- 
generate (single atomic specie) two-mode BEC problem, 
carry out the same particle-hole mapping, and divide the 
timescale by two. 

To conclude, the apparently complex dynamics of an 
ultracold Fermi gas coupled to a molecular BEC can be 
greatly simplified due to the inherent symmetry of the 
problem. We have obtained a system of dynamical equa- 
tions for N + 1 effective amplitudes enabling an improve- 
ment of three orders of magnitude in the total particle 
number of a realistic calculation. The structure of this 
system provides great insight into the origin of boson- 
like collective behavior. It demonstrates that dissocia- 
tion is enhanced as the number of particles whereas as- 
sociation is enhanced as the number of holes, due to con- 
structive interference of various pathways between states 
with different molecule number. Formal identity between 
fermion and boson dynamics was shown to exist when 
N = M and confirmed by numerical calculations. This 
equivalence can serve as an important tool in further 
studies of nearly degenerate fermion-boson systems. 
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